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REFLECTED BROWNIAN MOTION IN A CONVEX POLYHEDRAL CONE: 
TAIL ESTIMATES FOR THE STATIONARY DISTRIBUTION 

ANDREY SARANTSEV 


Abstract. Consider an multidimensional obliquely reflected Brownian motion in the positive orthant, 
or, more generally, in a convex polyhedral cone. We find sufficient conditions for existence of a stationary 
distribution and convergence to this distribution at the exponential rate, as time goes to infinity. We also 
prove that certain exponential moments for this distribution are finite, thus providing a tail estimate for 
this distribution. Finally, we apply these results to systems of rank-based competing Brownian particles. 


1. Introduction 

A multidimensional obliquely reflected Brownian motion in a convex polyhedron D C has been 
extensively studied in the past few decades. This is a stochastic process Z = {Z(t),t > 0) which takes 
values in D; in the interior of D, it behaves as a Brownian motion, and as it hits the boundary dD, it 
is reflected inside D, but not necessarily normally. For every face Di of the boundary dD^ there is a 
vector Ti on this face, pointing inside r*, which governs the reflection. If r, is the inward unit normal 
vector to D, then this reflection is normal] otherwise, it is oblique. Special care should be taken for the 
reflection at the intersection of two or more faces. A formal definition is given in Sect. 2. 

One particularly important case is the positive orthant D = R^, where R+ := [0,oo). The concept 
of a semimartingale reflected Brownian motion (SRBM) in the orthant was introduced in [18], as a 
diffusion limit for series of queues, when traffic intensity at each queue tends to one (the so-called 
heavy traffic limit). Later, it was applied in the theory of competing Brownian particles (systems of 
rank-based Brownian particles, when each particle has drift and diffusion coefficients depending on its 
current ranking relative to the other particles), see [21 [ 2 ]- The gap process (vector of gaps, or spacings, 
between adjacent particles) turns out to be an SRBM in the orthant. 

We refer the reader to the comprehensive survey [SB] about an SRBM in the orthant. Reflected 
Brownian motion in a convex polyhedron was introduced in [9] in a semimartingale form: semimartingale 
reflected Brownian motion^ or an SRBM. The paper [9] contains a sufficient condition for weak existence 
and uniqueness in law; it is stated in Sect. 2. 

In this paper, we assume that the condition mentioned above holds; then an SRBM in a convex 
polyhedron D exists and is unique in the weak sense, and versions of this SRBM starting from different 
points X G D form a Feller continuous strong Markov family. 

Of particular interest is a stationary distribution for an SRBM in a convex polyhedron D\ a prob¬ 
ability distribution vr on D such that if Z(0) ~ vr, then Z(t) ~ vr for all f > 0. This was a focus 
of extensive research throughout the last four decades. For the orthant (and, more generally, for a 
convex polyhedron), a necessary and sufficient condition for existence of a stationary distribution is not 
known. However, there are fairly general sufficient conditions and necessary conditions for the orthant, 
see For dimensions d = 2 and d = 3, a necessary and sufficient condition is actually known, see 

[27] and m Appendix A] for d = 2, and [1] for d = 3. For a convex polyhedron (more specifically, a 
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convex polyhedral cone), see ms for sufficient condition for existence of a stationary distribution. It 
was shown in [8] that if a stationary distribution exists, then it is unique. 

Exact form of this stationary distribution is known only in a few cases, the most important of which is 
the so-called skew-symmetry condition. Under this condition, the stationary distribution has a product- 
of-exponentials form, see [inii2n]. Other known cases (sums of products of exponentials) are studied 
in HD]. A necessary and sufficient condition for a probability distribution to be stationary is a certain 
integral equation, called the Basic Adjoint Relationship. However, it is not known how to solve this 
equation in the general case. 

We complement the results above by Ending some new conditions for existence of a stationary distri¬ 
bution for an SRBM in the orthant and in a convex polyhedral cone. To this end, we find a Lyapunov 
function: this is a function V : D ^ [1,00) such that for some constants k,b > 0 and a set C C D 
which is “small” in a certain sense (we specify later what this means; for now, it is sufficient to take a 
compact C) the process 

V{Z{t))-V{Z{0))- [ {-kV{Z{s)) + blc{Z{s)))ds 

Jo 

is a supermartingale. This is a more general definition than is usually used (with the generator of an 
SRBM). Under some additional technical conditions (so-called irreducibility and aperiodicity, more on 
this later), if such function V can be constructed, then there exists a unique stationary distribution 
TT, and the SRBM Z = {Z{t),t > 0) converges weakly to tt as t —>■ 00 ; moreover, the convergence is 
exponentially fast in t. There is an extensive literature on Lyapunov functions and convergence, see [21] 
for discrete-time Markov chains and [221 |23l ini HH HE] for continuous-time Markov processes. These 
methods were applied to an SRBM in the orthant in [T^ [6] and to an SRBM in a convex polyhedral 
cone in mi. However, in these articles they construct a Lyapunov function indirectly. In this article, 
we come up with an explicit formula: 

V{x) = U{x) := [x'QxY^^ x e D, 

where Q is a d x d symmetric matrix such that x'Qx > 0 for x E D \ {0}, A > 0 is a certain constant 
(to be determined later), ip : M_|_ —)■ M+ is a function with 


( 1 ) 


9?(s) 


0, s < Si; 

S, S > S2, 


for some 0 < Si < S 2 . 


We can also conclude that fj^V(x)7r(dx) < 00 . This explicit form of the function V allows us to find 
tail estimates for the stationary distribution tt. Let us also mention the paper im. which studies tail 
behavior of tt in case d = 2. 

The paper is organized as follows. In Sect. 2, we introduce all necessary concepts and definitions, 
explain the connection between Lyapunov functions, existence of a stationary distribution, and expo¬ 
nential convergence. In Sect. 3, we state the main result and compare it with already known conditions 
for existence of a stationary distribution; then, we prove this main result. Sect. 4 is devoted to systems 
of competing Brownian particles. 


1.1. Notation. We denote by R the k x fc-identity matrix, and by 1 the vector (1,..., 1)' (the dimension 
is clear from the context). For a vector x = (xi,... ,Xd)' G M'^, let ||x|| := {x\ + ... + be its 

Euclidean norm. The norm of a d x d-matrix A is defined by 


Pll 


max||Ax|| = max{\/A | A is an eigenvalue of A'A}. 

||a:||=l 


For any two vectors x,y E M'^, their dot product is denoted hj x ■ y = xiyi -f • • • -f Xdyd- As mentioned 
before, we compare vectors x and y componentwise: x < y if Xi < for all i = l,...,d; x < y 
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if Xi < Ui for all i = 1, ..., cf; similarly for x > y and x > y. This includes infinite-dimensional 
vectors from We compare matrices of the same size componentwise, too. For example, we write 
X > 0 for T G if Xj > 0 for f C = (cij)i<ij<d > 0 if Cij > 0 for all i, j. Fix 

d > 1, and let I C d} be a nonempty subset. Write its elements in the order of increase: 

I = {ii,... ,fm}, I < ii < i 2 <■■■< ira < d. For any x G W^, let [x]i := {xi^,.. .,Xi^y. For any 
d X d-matrix C = (cij)i<ij<rf, let [C]i := A one-dimensional Brownian motion with zero 

drift and unit diffusion, starting from 0, is called a standard Brownian motion. The symbol mes denotes 
the Lebesgue measure on We write / G C°°{D) for an inhnitely differentiable function / ; D —)■ M. 
Take a measurable space (X, z/). For any measurable function / ; X —)■ M, we denote (zz,/) := J^fdu. 
For a measurable function / : X ^ [1, cxd), dehne the norm ||z/||/ := sup \{u,g)\, where the supremum is 
taken over all measurable functions : X —)■ M such that \g{x)\ < f{x) for all t G X. For / = 1, this is 
the total variation norm: ||z/||tv- 


2. Deeinitions and Background 


2.1. Definition of an SRBM in a Convex Polyhedron. Fix the dimension d > 2, and the number 
m of edges. Take m unit vectors Ui,..., n^, and m real numbers &i,..., hm- Consider the following 
domain: 

D := {x E \ X ■ xii > bi, i = 1,..., m}. 

We assume that each face Di of the boundary dD: 

Di := {x E \ X ■ xii = bi, x ■ xXj > bj, j = 1,...,m, j 7 ^ f}, i = 1,... ,m, 

is (d — l)-dimensional, and the interior of D is nonempty. Then D is called a convex polyhedron. For 
each face Di, Uj is the inward unit normal vector to this face. Dehne the following d x m-matrix: 
N = [uil ... |nm]- Now, take a vector y eW^ and a positive dehnite symmetric d x d-matrix A. Consider 
also an m X d-matrix R = [ri |... |rm], with n* ■ r* = 1 for f = 1,..., m. We are going to dehne a process 
Z = {Z{t),t > 0) with values in D, which behaves as a d-dimensional Brownian motion with drift vector 
/i and covariance matrix A, so long as it is inside D] at each face it is rehected according to the 
vector ri. First, we dehne its deterministic version: a solution to the Skorohod problem. 


Definition 1. Take a continuous function X : M_|_ —)■ with X(0) E D. A solution to the Skorohod 
problem in D with reflection matrix R and driving function X is any continuous function Z : R+ —)■ 
D such that there exist m continuous functions £i,...,£m : R+ — )■ R+ which satisfy the following 
properties: 

(i) each is nondecreasing, £*(0) = 0, and can increase only when Z E Dy we can write the latter 
property as 

Ciit) = [ l{Z{s) E Di)dCi{s), t > 0; 
do 

(ii) for all t > 0 , we have: 


Z{t) = X{t) + RC{t), wheieCfl) = (£i(t),..., Cmit))'. 
The function £j is called the boundary term corresponding to the face D^. 


For the rest of the article, assume the usual setting: a hltered probability space (D, J-”, (d-))t>o, P) 
with the hltration satisfying the usual conditions. 


Definition 2. Fix z E D. Take an ((d-'t)i>o, P)-Brownian motion W = {W{t),t > 0) with drift vector 
y and covariance matrix A, starting from x. A continuous adapted D-valued process Z = {Z(t),t > 0), 
which is a solution to the Skorohod problem in D with reflection matrix R and driving function W, is 
called a semimartingale reflected Brownian motion (SRBM) in D, with drift vector y, covariance matrix 
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A, and reflection matrix R, starting from z. It is denoted by SRBM'^(Zi), R, /x, A). For the case D = 
we denote it simply by SRBM'^(R,/i, R). 

We shall present a snfficient condition for existence and nniqneness taken from j9]. First, let us 
introduce a concept concerning the geometry of the polyhedron D. 

Definition 3. For a nonempty subset I C { 1 ,..., m}, let Dj := riig/D,, and let D 0 := D. A nonempty 
subset I C { 1 ,..., m} is called maximal if Dj 7 ^ 0 and for I ^ J ^ {1,..., m} we have: Dj C Dj. 

Now, let us dehne certain useful classes of matrices. 

Definition 4. Take a d x d-matrix M = {mij)i<ij<d- R is called an S-matrix if for some u G m > 0 
we have: Mu > 0. It is called completely-S if for every nonempty I C {1 ,..., d} we have: [M]/ is an 
iS-matrix. It is called a ^-matrix if mij < 0 for i j. It is called a reflection nonsingular Ai-matrix if 
it is both a completely-iS and a ^-matrix with diagonal elements equal to one: rjj = 1, i = 1,..., d. It 
is called strictly copositive if x'Mx > 0 for x G \ {0}- It is called nonnegative if all its elements are 
nonnegative. 


A useful equivalent characterization of reflection nonsingular Ai-matrices is given in [2ll Lemma 2.3]. 
Now, let us hnally state the existence and uniqueness result, taken from [9|. 

Proposition 2.1. Assume that for every maximal set I C { 1 ,... ,m} the matrices [N'R]j and [RW]/ 
are S-matrices. Then for every z ^ D, there exists a weak version of an SRBM‘^(D, R, p. A), and it is 
unigue in law. Moreover, these processes for z E D form a Feller continuous strong Markov family. 

Remark 1. For a particular important case of the positive orthant: D = that is, when m = d, 
n* = Ci and 6 * = 0 for z = 1,..., d, we have: N = R, every nonempty subset / C {1,..., d} is maximal, 
and the condition from Proposition 12.11 is equivalent to the matrix R being completely S (because R is 
completely S if and only if R' is completely S). This turns out to be not just sufficient but a necessary 
condition, see [ 2 ^ . 


A sufficient condition for strong existence and pathwise uniqueness was found in [TS] for the orthant: 
R must be a reflection nonsingular-Af matrix. However, we shall not need strong existence and pathwise 
uniqueness in this paper. The generator of this process is given by 

i=l j=l 

with the domain R(A.) containing the following subset of functions: 

V{A) 0 {/ G C^{D) I r, ■ V/(x) 1,,^^ = 0, z = 1,..., m}. 


2.2. Recurrence of Continuous-Time Markov Processes. Let us remind the basic concepts of 
recurrence, irreducibility and aperiodicity for continuous-time Markov processes. This exposition is 
taken from [221 [231 [m E]. Take a locally compact separable metric space X and denote by 53 its 
Borel (T-held. Let 

(H,R, (R,) 

t>0) (A(f),f>0),(PR.gR 

be a time-homogeneous Markov family, where X{t) has continuous paths under each measure P^,. Denote 
by R*(x, A) = 'PxiXit) G A) the transition function, and by E^, the expectation operator corresponding 
to Rx- Denote by P^f and z/R^ the action of this transition semigroup on functions / : X —)■ M and 
Borel measures on X. Take a a-hnite reference measure v on X. The process X is called u-irreducible 
if for A G 53 we have: 


U(X(s))ds 


uo 


Z/(A) >0 Ea, 


> 0 for all X E X. 
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If such measure exists, then there is a maximal irreducibility measure V’ (every other irreducibility 
measure is absolutely continuous with respect to -0), which is unique up to equivalence of measures. 
A set A G with V’(^) > 0 is accessible. A nonempty (7 G is petite if there exists a probability 
distribution a on M_|_ and a nontrivial cr-hnite measure on 03 such that 

poo 

Wx eC, / P\x, ■)a{dt) > Ua{-). 

Jo 

Suppose that, in addition, this distribution a is concentrated at one point t > 0: a = 6t. Equivalently, 
there exists a t > 0 and a nontrivial a-hnite measure Ua on 03 such that 

poo 

Wx eC, / P\x, ■)a{dt) > Ua{-). 

Jo 

Then the set C is called small. The process is Harris recurrent if, for some a-hnite measure u, 


z/(A) > 0 ^ / lA(-A(s))ds = oo Pa, — a.s. for all x E X. 

Jo 

Harris recurrence implies z/-irreducibility. A Harris recurrent process possesses an invariant measure 
TT, which is unique up to multiplication by a constant. If tt is hnite, then it can be scaled to be a 
probability measure, and in this case the process is called positive Harris recurrent. An irreducible 
process is aperiodic if there exists an accessible petite set C and T > 0 such that for all x G C and 
t >T, we have: P*(x, C) > 0. 


Definition 5. The process X is called V-uniformly ergodic for a function H : X —)■ [l,cxo) if it has a 
unique stationary distribution tt, and there exists constants K, > 0 such that for all x G X and t >0 
we have: 


P*(x,-) - vr(-)|| < KV{x)e-^. 


Now, let us state a few auxiliary statements. The next proposition was proved in [211 Chapter 6] for 
discrete-time processes, but the proof is readily transferred to continuous-time setting. 


Proposition 2.2. For a Feller continuous strong Markov family, every compact set is petite. 

Lemma 2.3. Take a Feller continuous strong Markov family. Assume t/j is a reference measure such 
that there exists a compact set C with > 0. If P^{x, A) > 0 for all t > 0, x E X and A G 53 such 
that '0(A) > 0, then the process is tlr-irreducible and aperiodic. 

Proof. Irreducibility follows from the dehnition. For aperiodicity, we can take the compact set C, 
because it is petite by Proposition 12.21 If ip' is a maximal irreducibility measure, then friC) > 0 and 
0 -C Tp', and so Tp'iC) > 0. The rest is trivial. □ 

Finally, the following statement was proved in [HI Lemma 3.4]. 

Proposition 2.4. For an SRBM‘^(P, P, p. A) under the conditions of Proposition I^TTl for every t > 0, 
X E D, and A D with mes(A) > 0 we have: P*(x, A) > 0. 

Remark 2. Combining Lemma [2731 and Proposition 12.41 we get that an SRBM'^(P, R, /x. A) is irreducible 
and aperiodic. 

2.3. Lyapunov Functions and Exponential Convergence. There is a vast literature (some of these 
were mentioned in Sect. 1) on connection between Lyapunov functions for Markov processes and their 
convergence to the stationary distribution. However, for the purposes of this article, we need to state 
the result is a slightly different form. First, let us dehne the concept of a Lyapunov function. 
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Definition 6. Take a continuous function V ■. X ^ [l,C)o). Suppose there exists a closed petite set 
C C X and constants fc, 6 > 0 such that the process 

(2) V{X(t)) - l"(X(0)) - f [-kV{X{s)) + blc{X{s))] ds 

Jo 

is an ((J^t)i>o, Pa;)-superniartingale for all x G X. If, in addition, sup^-P < oo, then V is called a 
Lyapunov function for the process X. 

Remark 3. Equivalently, we can request that the process in ([2]) is a local supermartingale. This is 
equivalent to it being a supermartingale, because this process is bounded from below by —V{x) — bT 
on any time interval [0,T] under the measure Pa,. (Every local supermartingale which is bounded from 
below is a true supermartingale; this follows from a trivial application of Eaton’s lemma.) 

This definition is taken from m Section 3] with minor adjustments, with (^(s) = ks in the notation 
of [m. This is a slightly more general definition than is often stated in the literature; a more customary 
one invloves the generator of the Markov family. First, let us state an auxiliary lemma. 

Lemma 2.5. For some constant cq > 0, we have: P^U{x) < c^U{x), for all x & X and s G [0,1]. 

Proof. Because U and V are equivalent in the sense of dl]), it suffices to prove the statement of Lemma 1231 
for V instead of U. But this follows from the fact that the process ([2]) is a super mart ingale. Indeed, 
take Ej; in ([ 2 ]) and get: 

PW{x)-V{x) + k [ P^V{x)-b[ P*(x,C')ds < 0. 

Jo Jo 

Therefore, P^V{x) < V{x) + bt. But V{x) > 1, so for t G [0,1], we get: P^V{x) < (1 + b)V{x). This 
completes the proof of Lemma 12.51 □ 

Next, we present the main result for this subsection. 


Theorem 2.6. Assume there exists a Lyapunov function V, and the process is irreducible and aperiodic. 
Then there exists a unique stationary distribution vr, the process is V-uniformly ergodic, and we have 
the following estimate: 


( 3 ) 


(vr, P) = / I/(x) 7 r(dx) < oo. 


Proof. Existence and uniqueness of tt together with ([3]) follows from m Proposition 3.1]. If the process 
is irreducible, then the skeleton chain {X(n))nez+ is irreducible. Apply |TTl Theorem 3.3] to the case 
(p(x) = fcx, we get that for any to > 0, there exists a function V : X ^ [k, oo), an accessible petite set 
C for the skeleton chain {X{n))nei,+ and a constant 6 > 0 such that sup (5 P < oo. 


0 < Cl < 


V{x) 

V{x) 


< C2 < oo. 


xgX, and < (1 - fc)P + 61 ^. 


Taking U := P/fc : X —)■ [l,oo), we get: there exists A := 1 — /c < 1 and b' = b/k > 0 such that 
P^U < -XU + b'lc, and 


( 4 ) 


0 < C3 < 


P(x) 

V{x) 


< C4 < oo. 


X G X. 


It follows from Proposition 12.41 that the skeleton chain (X(?7,))„>o is irreducible and aperiodic. By [211 
Theorem 5.5.7], the petite set C is small for this skeleton chain. Since this chain is irreducible and 
aperiodic, by [121 Theorem 2.1(c)] for some constants C5 > 0 and p G (0,1), we have: 


(5) 


P"(x,-)-vr(-)|k<C5P(x)pT 
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Next, we follow the proof of [T^ Theorem 5.2], Every f > 0 can be represented as t = n+s, where n G Z_|_, 
s G [0,1). Since vr is stationary, we have: vrP® = vr. Therefore, for any measnrable g : X ^ [1, oo) with 
|c/(t)| < U{x), 

P*g{x) - {7r,g) = P^P^g{x) - {7r,P^g). 

Bnt from Lemma [2.51 we have: 


\P"g{z)\ < P^u{z) < ceU{z), zeX. 
From ([5]), becanse \g{x)\ < U{x) for x E X, we get: 

\P^P^g{x) - {n,P^g)\ < c,c^U{x)p^. 


Since n < f — 1, 

\P^g{x) - {'x,g) \ < c^c&p~^U{x)pK 
This proves that, for cy := CsCep”^ 

\\P\x, •) - 7 r (-)||;7 < C 7 U{x)e~’^, K := -Inp. 

This is t/-nniform ergodicity. Since the fnnctions P and V are eqnivalent in the sense of (jl]), this also 
means E-nniform ergodicity. □ 


3. Main Results 


3.1. Statement of the General Result. Consider now a special type of a convex polyhedron, namely 
a convex polyhedral cone: P = {x G | Nx > 0}, where N is a m x d-matrix, constrncted in Sect. 

2.2. This hts into the general framework of Dehnition [2|, if we let 6 i = ... = 6 ^ = 0. What follows is 
the main resnlt of the paper. 


Theorem 3.1. Suppose that conditions of Proposition 12.11 hold, 
nonsingular d x d-matrix Q such that: 

(i) x'Qx > 0 for X E D \ {0}; 

(a) (R'Qx)i < 0 for X E Dj, for each ? = l,...,m; 
(zzz)x'Qp<0forxED\{0}. 

Take a C°° function (p : R+ —)■ R+ defined in (PP). Denote 


( 6 ) 


A : = 


2 min 

3;GD\{0} 


\Qp ■ x\U{x) 
x'QAQx 


Assume there exists a symmetric 


Then for A G (0, A), the function 

(7) Vx{x) = U{x) := [x'Qxf^ , 

is a Lyapunov function for the SRBM'^(Z1, P, /i, A). Therefore, the SRBM'^(Z1, P, p. A) has a unique 
stationary distribution n, which satisfies 


( 8 ) 

and is V\-uniformly ergodic. 


(71,14)= / 14(T)7r(da:) < oo, 
J D 


Remark 4. The qnantity A is strictly positive. Indeed, the matrix A is positive dehnite, and Q is 
nonsingnlar; so for a: 7 ^ 0 we have: Qx 7 ^ 0 and x'QAQx = {QxfA^Qx) > 0. Also, Qp ■ x = x'Qp < 0, 
and U{x) > 0 for a: G D\ {0}. Therefore, the fraction is positive for each x E D \ {0}. Since 
this fraction is homogeneons (invariant nnder scaling), we can take the minimnm on the compact set 
{x E D \ 1411 = 1}. The rest is trivial. 
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The estimate ([8]) implies that some exponential moments of tt are hnite. Namely, let 


(9) 


K := min Uix). 

xeD 

\\x\\=l 


This quantity is strictly positive, because U{x) >0 on the compact set {x G -D | ||x|| = 1}. Therefore, 
for large enough ||x|| we have: 

Vx{x) > 

and 

f < oo for p E {0,AK). 


' D 


From here, we get: for every a > 0, 


7r{x G D I ||x|| > a} < C{p)e for p G (0, AiF). 

Let us compare this result with DEI. where a more general case is considered (drift vector and 
covariance matrix depend on the state). There, a sufficient condition for Id-uniform ergodicity is: 

(i) that Skorohod problem in D has a unique solution for every driving function and is Lipschitz 
continuously dependent on this function, in the metric of C'([0, T], M'^) for every T > 0; 

(ii) there exists a vector h G M™, b > 0, such that Rb = —p. 

Condition (i) is stronger than the one from Proposition 12.11 However, we were not able to come up 
with an example when conditions of Theorem 13.11 hold, but condition (i) does not hold. Some sufficient 
conditions for (i) to hold are known from [13]. However, condition (ii) is much simpler than (i) - (hi) 
from Theorem |3.11 The results from DE] also construct a Lyapunov function indirectly, without giving 
an explicit formula. This does not allow to construct explicit tails estimates, as in 13.11 


3.2. Applications to the Case of the Positive Orthant. Now, let D = that is, m = d and 

N = Id- We have the following immediate corollary of Theorem 13.11 


Corollary 3.2. Assume R is a completely S matrix. Suppose there exists a strictly copositive nonsin¬ 
gular d X d-matrix Q such that QR is a Z -matrix, and Qp < 0. Then an SRBM'^(i?, p, A) has a unique 
stationary distribution n and is V\-uniformly ergodic for X G (0, A), while tt satisfies ([H]). Here, Vx is 
defined in ([7]), and A is defined in (|2]). 

Proof. Condition (i) of Theorem 13.11 follows from the definition of copositivity. Condition (ii) follows 
from the assumption that QR is a 2^-matrix, because then for z E Di we have: z > 0, but Zi = 0, and 
so 

d 

{QRz)i = '^{QR)ijZj = '^{QR)ijZj < 0. 
i=i 

Condition (iii) follows from Qp <0. □ 

A particular example of this is as follows. 


Corollary 3.3. Assume R is a d x d-refilection nonsingular AA-matrix, and there exists a diagonal 
matrix C = diag(ci, ... ,Cd) with ci,..., q > 0 such that R = RC is symmetric. If R~^p < 0, then the 
process SRBM'^(i?, p, A) has a unique stationary distribution tt, and is Vx-uniformly ergodic with 


\4(a;) = U{x):= 

for X G (0, A), where the function ip is defined in ([T]), and 

A := 2 


■ /—1 1 ^/^ 
X R X 


mm 


R ^p ■ x\U{x) 


9 ^GK^\{o} x'R ^AR \ 
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In addition, (vr, V\) < oo for A G (0, A). 

Proof. Just take Q = R ^ = C~^R~^ in Corollary 13.21 Let us show that the matrix Q is strictly 
copositive. From [2H Lemma 2.3], R~^ is a nonnegative matrix with strictly positive elements on the 
main diagonal. Since C~^ is a diagonal matrix with strictly positive elements on the diagonal, the matrix 
R is also a nonnegative matrix with strictly positive elements on the main diagonal. Therefore, for 
X G X 7^ 0 we have: x'R x > 0. Now, from < 0 it follows that R /i < 0. The rest is 

trivial. □ 


Example 1. However, Corollary 13. 2l can be applied not only to the case when i? is a reflection nonsingular 
AJ-matrix. Indeed, let d = 2 and 


R 


1 0.5 

0.5 1 


/i 



Then i? is a completely S matrix. Take the matrix 


Q 


1 - 0.6 

- 0.6 1 


then QR 


0.7 -0.1 

-0.1 0.7 


is a 2^-matrix, and Q/a < 0. However, R is not a reflection nonsingular AJ-matrix. 


It is instructive to compare these results with already known ones. It turns out that the only new 
statement in Corollary 13.31 is the tail estimate (tt, Va) < cxo for an explicitly constructed function V\. 
Existence (and uniqueness) of a stationary distribution and H-uniform ergodicity for some function 
V : —>■ [l,cx3) are already known from [15], [6] (only there is no simple formula for the Lyapunov 

function V in these papers: it is simply known that V{x) > for some 01,02 > 0). The paper 

[15] states the fluid path condition, which is sufficient for H-uniform ergodicitiy: for every x G M^, any 
solution of the Skorohod problem in the orthant with reflection matrix R and driving function x + pit 
must tend to zero as t —)■ 00 . This turns out to be a necessary and sufficient condition for the case 
d = 3. In the case d = 2, another necessary and sufficient condition is found: R must be nonsingular 
and < 0; see [27] and m Appendix A]. In fact, the following condition is necessary for existence 

of a stationary distribution: R is nonsingular and R~^fi < 0, see [H Appendix Cj. For d = 3, the fluid 
path condition is weaker than this necessary condition, see [1]. 

It is not known for d > 4 whether the fluid path condition is necessary. Therefore, Corollary 13.21 
might contain results which are new compared to the fluid path condition. However, we do not know 
any counterexamples to fluid path condition (that is, cases when it is false, but the stationary distribution 
exists). This is a matter for future research. 

3.3. Proof of Theorem 13.11 Recall from Remark [2] that an SRBM'^(Z1, R, p. A) is irreducible and 
aperiodic. The rest of the proof will be devoted to proving that the function (I7|) is indeed a Lyapunov 
function in the sense of Dehnition [6l Apply the Ito-Tanaka formula to 

Z{t) = W{t) + RL{t), t>0, 

where Z = {Z{t),t > 0) is an SRBM'^(iA, R, p. A), W = {W{t),t > 0) is the driving Brownian motion 
for Z, and L = (Li,... ,Lm)' is the vector of boundary terms. Because of (II]), for x G R such that 
||x|| > S 2 , we have: 

Vy,{x) = 

First, let us calculate the hrst- and second-order partial derivatives of U on this set. Since Q is symmetric 
and x'Qx > 0 for x G R such that ||x|| > S 2 , we have: 

d{x'Qx) 
dxi 


2(^Qxf, % 1,..., d. 
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Therefore, 

dU{x) 

1 dU{x) 

{Qx)i , 


dxi 

2U (x) dxi 

~ U{x) ’ 

Now, 

S^uix) A§Au{x) 

- (0-).^ 

qijU{x) - {Qx)i 


U{x) 


dxidxj 


U^{x) 

qijU^{x) - {Qx)i{Qx) 


U^x) 

{qij{x'Qx) - (Qxx'Q)ij). 


U^{x) U^{x) 

As ||a;|| oo, these second-order derivatives tend to zero, because U{x) > K\\x\\ for x E D. 
us calculate the first- and second-order partial derivatives for Vx'. 

dVx{x) dU ^{Qx)i 


and 


d‘^Vx{x) _ ^dVx{x) {Qx)i , ^^^{Qx)i{Qx)j , ,^-,d‘^U{x] 


dxidxj 


= \- 


dxj U (x) 


+ AVx(a:) 


dxidxj 


= yvx{x)- 


U{x) U{x) 


+ XVx{x) 


dxidxj 


_ ^2y (AQ^^'Q)ij , ,^U{X) 

- A ^a(x) + ^Vx{x)^^ 


x'Qx 


Since (Zj, Zj)t = aijdt for i, j = 1 ,..., d, and Z{t) = Id^(t) -|- RL(t), we have: 

dv(z W) = i E E +1 


2 dxidxj 

1=1 j=i * J 

^EE“« (\^Vx{Z{t)) 


2 = 1 


dxi 


i=l j=l 
d 


{QZ{t)Z{t)'Q) 


V 




{Qzm 

U{Z{t)) 


/i. 


^-Vx{Z{t))dt + J2^ 


Z{t)'QZ{t) 

QZ{t) 

j 


+ \Vx{Z{t)) 


^ U{Z{t)) 


d^U{Z{t)) \ 

dxidxj J 

Vx{Z{t))r^dL^{t) 


dt 


+ AV(Z(«)) E (»■.(«) - 

= V(Z(i))/?i(Z(i))d« + dM(«) + dL(«), 

where for a; G -D \ {0} we let I3x{x) := a(a;)A^ — h{x), where 

{Qxx'Q)ij tr{AQxx'Q) tT{x'QAQx) x'QAQx 


d d 


«(^) = EE' 


Hj 


i=l j=l 


x'Qx 


x'Qx 


and 


_6(^) e{x) + where 9(x) ^ E E “« 

^ ^ i=l j=l * •? 


x'Qx x'Qx 

d^U{x) 


d d 


and, in addition. 


M{t):=\ / Vx{Z{t))Y,A 


mm 

U{Z{t)) 


d{Wi{t) - Hit ), 


Now, let 
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m Evilly'' a(Z(*)) • r,ALi{t). 

Lemma 3.4. The process L = {L{t),t > 0) is nonincreasing a.s. 

Lemma 3.5. For A < A, there exist r(A),/c(A) > 0 such that for x E D, ||a;|| > r(A), we have: 
Pxix) < -k{X). 

Assuming we proved these two lemmata, let us complete the proof of Theorem 13.11 Fix A G (0,A). 
Take a compact set C = {x G -D \ ||a;|| < 'r(A)}, with r(A) from Lemma [3.51 By Proposition 12.21 this 
set is petite. The process 

Vx{Z{t)) - Vx{Z{0)) - [ Vx{Z{s))/3x{Z{s))ds, t > 0, 

Jo 

is a local supermartingale, because M = {M{t),t > 0) is a local martingale and by Lemma [3.41 Now, 

/3x{x)Vx{x) < -k{X)Vx{x)lD\c{x) +bxlc{x), 

for bx ■= maxa,g (7 [(3x{,x)Vx{x)]. This maximum is well dehned, because (JxVx is a continuous function, 
and C is a compact set. The rest of the proof is trivial. 

Proof of Lemma lATl We can write L(t) as 

dpt) = 

But each Lj can grow only when Z G Dj, and then {R'QZ{t))j = QZ{t) - rj < 0. It suffices to note that 
Vx{Z{t)) > 0 and U{Z{t)) > 0. 

Proof of Lemma [331 For each x E D \ {0} we have: if b{x) > 0, then 

A < A(a:) := \ | /Jx{x) < 0. 
a[x) 

Note that d{x) —)■ 0 as ||a;|| —)■ cx3. From this and conditions (i), (ii) and (iii) of Theorem 13.11 it is 
straightforward to see that 

lim A(a;) = A. 

||a;||^oo 

x£D 

Also, there exist ro,co > 0 such that for x E D, ||x|| > vq we have: a{x),b{x) > cq. Now, £x A G (0,A). 
Then there exists 5 > 0 such that S < X < A —26, and there exists r(A) such that for x E D, ||x|| > r(A) 
we have: A(a;) > A — 6. Without loss of generality, we assume r(A) > vq. Now, for such x we have: 

—/3x{x) = —a{x)X‘^ + b{x)X = a(a;)A(—A+A(a;)) > Co5((A —5) — (A —25)) > CqS'^. 

This completes the proof of Lemma 13.51 

4. Systems of Competing Brownian Particles 

4.1. Classical Systems: Definitions and Background. In this subsection, we use definitions from 
[3]. Assume the usual setting: a hltered probability space (fl, 3, {J^t)t>o, P) with the hltration satisfying 
the usual conditions. Let N > 2 (the number of particles). Fix parameters gi,...,g]^ E M and 
(Ti,..., (Ttv > 0. We wish to dehne a system of N Brownian particles in which the kth smallest particle 
moves a Brownian motion with drift g^ and diffusion cr^. Applications include: (i) mathematical finance, 
namely modeling the real-world feature of stocks with smaller capitalizations having larger growth rates 
and larger volatilities; it suffcies to take decreasing sequences (gk) and (u^); (ii) diffusion limits of a 
certain type of exclusion processes, namely asymmetrically colliding random walks, see 1211. 
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Definition 7. Take i.i.d. standard (J^t)t>o-Brownian motions Wi,..., Wj^. For a continuous M^-valued 
process X = {X(t), t > 0), X{t) = (Xi(t),... ,Xjv(t))', let us define p^, t > 0, the ranking permutation 
for the vector X{t): this is a permutation on {1, ..., iV}, such that: 

(i) Xp,p)(t) < Xp,(^)(t) for 1 < i < j < TV; 

(ii) if 1 < i < j < X and Xp,p)(f) = Xp,(j)(f), then pt{i) < Pt{j)- 
Suppose the process X satisfies the following SDE: 

N 

(10) dXi{t) = '^l{pt{k) =i)[gkdt + akdWi{t)], i = l,...,N. 

k=l 

Then this process X is called a classical system of N competing Brownian particles with drift coefficients 
gi, ... ,gN and diffusion coefficients af, ... ,ajf. For i = 1, ... ,N, the component Xj = (Xj(t), t > 0) is 
called the ith named particle. For k = 1, ..., N, the process 


Bfc — (Xfc(t), t > 0), Yk{t) Xpj(fc)(t) = X(fc)(t), 

is called the kth ranked particle. They satisfy Yi(t) < Y 2 {t) < ... < F/v(f), t > 0. If pt{k) = i, then we 
say that the particle Xft) = Yk{t) at time t has name i and rank k. 


Weak existence and uniqueness in law was established in |3]. Consider the gap process: an 
valued process defined by 


Z = iZit),t>0), = Zkit)=Yk+iit)-Ykit). 


It was shown in |3] that this is an SRBM^ ^{R,pi,A) in the orthant S = ^ with parameters 


( 11 ) 


( 12 ) 



■ 1 

-1/2 0 0 

... 0 O' 



-1/2 

1 -1/2 0 

0 0 



0 

-1/2 1 0 

0 0 


R = 






0 

0 0 0 

1 -1/2 



0 

0 0 0 

... -1/2 1 



T = {92 

~ 9ly 93 ~ fi'45 ■ ■ 

) 9n — 9n-i)' , 


'af 

-(^2 

0 0 

0 

0 

-0-2 cr 


-a| 0 

0 

0 

0 


al + aj -aj 

0 

0 

0 

0 

0 0 

2 2 
• • '^N-2 + '^N-1 


0 

0 

0 0 

—CTn-I 

+ a%_ 


4.2. Main Results. In this subsection, we present results about the gap process: existence of a station¬ 
ary distribution, Lyapunov functions and tail estimates. Let gg. := {gi + ... + gk) //c for fc = 1,..., X. 

Proposition 4.1. The gap process has a stationary distribution if and only if 

(13) 9k>9N for/c= 1,...,X-1. 

In this case, it is V-uniformly ergodic with a certain function V : —)■ [1, oo). 

Proof. This result was already proved in [3l IH [25] , but for the sake of completeness we present a sketch 
of proof. The matrix i? is a reflection nonsingular Al-matrix, and 

(14) - R”V = 2 (5-1 - 9n: 9i + 92- 2g^,..., gi + ... + g^-i - (X - . 
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Define the quantities 

( 15 ) h = gi + g2 + ... + gi-igj^, i = -1. 

Then we can rewrite flT4)) as 

-R~^H = 2b, b = {bi,... ,bN-iy. 

Therefore, the gap process has a stationary distribution if and only if each component of this vector 
is strictly positive, which is equivalent to the condition flT^ . In this case, the fluid path condition 
holds by [7], and so by [6] the gap process is D-uniformly ergodic for a certain Lyapunov function 
V ^ [l,oo). □ 

From Corollary 13.31 we get a concrete Lyapunov function V, namely: 

Vy^{x) = U{x) := [x'R-^x]^'^ , 

where ip is dehned in ([T]). We use the fact that the matrix R is symmetric, so in the notation of 
Corollary 13.31 we have: C = In-i and R = R. Here, we must have A < A, where 

A o • \R~^fi ■ x\U{x) 

x'R ^AR ^x 

Let us try to estimate the tail of the stationary distribution vr. 


Theorem 4.2. Using the definition of bi,... ,6Ar_i from flTSD . we have: 




gPlhll 7 ^(da;) < oo for p e (0,po), Po ■= 


2 min(6i,.. .,bN-i) 




Pll 


N- 


Proof. From the results of Theorem 13.11 we have: 

f < oo for pG (0,iCA), 


where K is dehned in ([9]) (in this notation, Q = R ^). Now, let us estimate KA from below. Dehne 
S := {x e I xi + ... + Xn-1 = !}• 

Lemma 4.3. (i) The norm of the matrix R~^ is equal to 

(16) ||i?-i|| = Ah^ = (^1-cos^) ; 

(a) U{x) >1 for X G S; 

(Hi) \R~^p • tI > 2 min(&i,..., bp^^i) for x G S; 

(^v) x'R-^AR-^x < \\R-ymA\\ forxe S. 


Suppose we proved LemmaFrom part (ii) we get: iF > 1. Using (ii)-(iv), we obtain: 

^ 2 min(&i,...,&jv_i) 

\\r-T\\M ■ 


Finally, using (i), we get: 

But 

Therefore, 


A > 4 1 


vr \2 min(6i,..., fejv-i) 
cos — — 


NJ 


Pll 


1 - 


71 

cos-> 


1 

2 



TT^ 1 

YJp' 


A > 


2 min(6i,...,6jv-i) ^,-2 

.2 ||H|| 
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The rest of the proof is trivial. 

Proof of Lemma IT731 (i) The eigenvalues of R are given by (see, e.g., [?]) 

/CTT 

Afc = l-cos—, k = 1,...,N -1. 

The eigenvalues of R~^ are fc = 1,..., iV — 1. The matrix R~^ is symmetric, so its norm is equal 
to the absolute value of its maximal eigenvalue. Therefore, we get (fT6D . 

(ii) The matrix R~^ is symmetric and positive dehnite. Solving the optimization problem x'R~^x —)■ 
min, a; • 1 = 1, we get: the minimum is l'i?l, which is equal to the sum of all elements of R, which, in 
turn, equals 1. 

(hi) Follows from the fact that R~^fi < 0 and flTTl) . 

(iv) Follows from the multiplicative property of the Euclidean norm, and from the fact that for x G S 
we have: ||x||^ = x^ + ■ ■ ■ + x%_i < (xi +-h x%_iY = 1. □ 

4.3. Asymmetric Collisions. One can generalize the classical system of competing Brownian particles 
from De£nition[7]in many ways. Let us describe one of these generalizations. Consider a classical system 
of competing Brownian particles, as in Definition [71 For fc = 1,..., iV — 1, let L(^k,k+i) = (7^(fc,A;+i)(^)) ^ — 
0) be the semimartingale local time process at zero of the process = Y^+i — Yk- We shall call this 
the collision local time of the particles Yk and W+i- For notational convenience, let L(o,i)(t) = 0 and 
L{N,N+i){t) = 0. Let 

Bk{t) = j l(ps(/c) = i)dlFi(s), /c = l,...,A, t > 0. 

i=i 

It can be checked that {Bk, Bi)t = so Bi,..., Bj\f are i.i.d. standard Brownian motions. As shown 
in pig, the ranked particles Yi,, Y/v have the following dynamics: 

Yk{t) = Ffc(O) + Qkt + akBkif) — -L(fc^fc_|_i)(t) + A: = 1,..., A. 

The collision local time L(^k,k+i) has a physical meaning of the push exerted when the particles Yk and 
Yk+i collide, which is needed to keep the particle Yk+i above the particle Yk. Note that the coefficients 
at the local time terms are ±1/2. This means that the collision local time L(^k,k+i) is split evenly between 
the two colliding particles: the lower-ranked particle Yk receives one half of this local time, which pushes 
it down, and the higher-ranked particle Y^+i receives the other one half of this local time, which pushes 
it up. In the paper [20], they considered systems of Brownian particles when this collision local time is 
split unevenly: the part qfj^]^L{k^k+i){i) goes to the upper particle Ffc+i, and the part L(^k,k+i){t) goes 
to the lower particle Yk. Let us give a formal definition. 

Definition 8. Fix N > 2, the number of particles. Take drift and diffusion coefficients gi,..., E M., 
(Ti,..., (Tat > 0, and, in addition, take parameters of collision 

Aw e (0)1)) ?a!+i + C = 1) A: = 1,...,A-1. 

Consider a continuous adapted M-^-valued process Y = (Y(t) = (Yi(t),... ^Yjsiit))'R > 0). Take other 
A — 1 continuous adapted real-valued nondecreasing processes 

B{k^k+i) (7^(fc,fc+i)(^)) ^ — 0), k 1,..., A 1, 

with L(fc fc_|_i)(0) = 0, which can increase only when Yk+i = Yk'. 

liYk+iit) > Yfc(A))dL(fc,fc+i)(A) =0, A: = 1,..., A - 1. 
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Let L(o,i)(t) = 0 and L( 7 v,Ar+i)(t) = 0. Assume that 

(17) Yk{t) = Yfc(O) + gkt + (TkBkit) - q'^ L(^k,k+i){t) + L(fc_i,fc)(t), k = l,...,N. 

Then the process Y is called the system of competing Brownian particles with asymmetric collisions. 
The gap process is dehned similarly to the case of a classical system. 

Strong existence and pathwise uniqueness for these systems are shown in [2DI Section 2.1], When 
gf = = ... = 1/2, we are back in the case of symmetric collisions. 


Remark 5. For systems of competing Brownian particles with asymmetric collisions, we dehned only 
ranked particles Yi, • • •, Y^- It is, however, possible to dehne named particles Xi,..., Xat for the case of 
asymmetric collisions. This is done in [201 Section 2.4]. The construction works up to the hrst moment 
of a triple collision. A necessary and sufficient condition for a.s. absence of triple collisions is given in 
[24] . We will not make use of this construction in our article, instead working with ranked particles. 


It was shown in [20] that the gap process for systems with asymmetric collisions, much like for the 
classical case, is an SRBM. Namely, it is an SRBM^“^(i?,/i. A), where p, and A are given by ffT2l) 
and ffT3]) . and the rehection matrix R is given by 


(18) 


■ 1 

-92 

0 

0 

0 

0 ■ 

-4 

1 


0 

0 

0 

0 

-4 

1 

-c 

0 

0 

0 

0 

0 

0 

1 

T 

.. 1 

1 

0 

0 

0 

0 

1 

1 

1 


This matrix is also a rehection nonsingular A4-matrix. Therefore, there exists a stationary distribution 
for this SRBM if and only if R~^p < 0. In this case, we can apply the results of [6] again and conclude 
that the gap process is R-uniformly ergodic with a certain Lyapunov function V : —)■ [1,cxd). 

Corollary 13.31 allows us to hnd an explicit Lyapunov function and provide explicit tail estimates. A 
remark is in order: the matrix R in flTSD in general is not symmetric, as opposed to the matrix R 
in (ITT]) . But for the following (X — 1) x (X — 1) diagonal matrix C, the matrix R = RC is diagonal: 


C 


diag 


^ 44 

’ (l2 ’ <l2 % ’' 


■■■ 4-1 \ 

(l2% ■■■(& -I ) ' 
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